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Abstract 

Given an equilateral triangle with a the square of its side length and a point in its plane 
with b, c, d the squares of the distances from the point to the vertices of the triangle, it can 
be computed that a, b, c, d satisfy 3{a^ + b'^ + + d'^) = {a + b + c + dY. This paper derives 
properties of quadruples of nonnegative integers {a, b, c, d), called triangle quadruples, satis- 
fying this equation. It is easy to verify that the operation generating {a, b, c, a + b + c — d) 
from {a,b,c,d) preserves this feature and that it and analogous ones for the other elements 
can be represented by four matrices. We examine in detail the triangle group, the group with 
these operations as generators, and completely classify the orbits of quadruples with respect to 
the triangle group action. We also compute the number of triangle quadruples generated after 
a certain number of operations and approximate the number of quadruples bounded by charac- 
teristics such as the maximal element. Finally, we prove that the triangle group is a hyperbolic 
Coxeter group and derive information about the elements of triangle quadruples by invoking 
Lie groups. We also generalize the problem to higher dimensions. 

1 Introduction 

The study of Apollonian circle packings is an elegant geometric subject that generates deep ques- 
tions in the theory of hyperbolic Coxeter groups, Lie theory, and analytic number theory, and many 
interesting problems about these packings can be effectively solved with state-of-the-art techniques 
from these fields. For example, [|2l|, [[31, and other papers in the same series by these authors de- 
rive a large number of properties of the Apollonian group and Apollonian quadruples from group 
theory, number theory, and geometry. The case of Apollonian equilateral triangles is a similar but 
less well-studied topic, and the goal of the project was to answer natural questions about this case 
with advanced techniques in group theory and number theory. 

Formally, an Apollonian circle packing is a fractal generated from triples of mutually tangent 
circles. In such a packing, for any four mutually tangent circles, their curvatures, a, b, c, d, satisfy 
Descartes's equation, 2{a^ + + c'^ + d^) = (a + b + c + d)^. Note that for fixed a, b, c, exactly two 
values of d satisfy the equation, generating two possible configurations, as illustrated in Figure [T] 
Given the configuration in Figure [Taj additional circles can be inscribed in each lune successively, 
a process illustrated in Figure |2j In addition, it is clear that given four mutually tangent circles with 
curvatures a, b, c, d, the curvature of the new circle inscribed in the lune bounded by the circles 
with curvatures a, b, c is 2(a + b + c) —d, which is easily derived from Descartes's equation. 
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(a) If the fourth circle is internally tangent to (b) If the fourth circle is externally tangent to 
the other three, then its curvature d is negative. the other three, then its curvature d' is positive. 



Figure 1 : There are two ways in which a fourth circle can be packed, and the two curvatures are related according to 

d^-d' =2(fl + fe + c). 




Define an Apollonian quadruple be a quadruple of nonnegative integers satisfying Descartes's 
equation. Then the inscribing operation described above can be expressed as the transformation 
(a, b, c, d) — )■ (a, b, c, 2(a -\-b + c)—d). Now, consider the matrices 

/I \ 
10 
10' 
\2 2 2 -1/ 

It is clear that for an Apollonian quadruple v = (a, b, c, dY , we have 54V = (a, Z?, c, 2(a + + 
c) — dY , and analogous relations hold for 5i, 5*2, S^,. Therefore, 5i, S2, ^3, ^4 correspond to 
transformations on a, b, c, d, respectively. Define the Apollonian group to be the group with 5i, 
52, ^3, 54 as generators. Previous papers have proved that the Apollonian group is a Coxeter group. 

Now, consider an equilateral triangle T and a point P in its plane. Let a denote the square of 
the side length of T, and let b, c, and d denote squares of the distances from P to the vertices of 
T. With simple trigonometry, we can derive the relation 3{a^ + b'^ + c'^ + d^) = (a + b + c + dY, 
which is remarkably similar to Descartes's equation. For the relation for the general case of an 
n-dimensional simplex, see [5J. 

In the Apollonian triangle case, as in the Apollonian circle case, there exist four matrices 
corresponding to reflection operations (instead of inscribing operations), as illustrated in Figure |3} 
with which a triangle group can be defined. Because of the great similarities between this case and 
the Apollonian case, we are motivated to verify whether many of the properties in the Apollonian 
case also hold for the triangle case, and where they do not, to determine whether the conclusions 
can be modified. We also analyze properties of the triangle group and triangle quadruples that have 
not been studied in the Apollonian case. 

In Section|2} we define the terms relevant to the triangle group and quadruples. We also discuss 
the reflection operation in greater detail. In Section [3} we introduce the algorithm of reducing 
triangle quadruples, which was also extensively analyzed for Apollonian packings. In Section]?} 
we prove that the triangle group is a Coxeter group. In Section[5| we compute several properties of 
triangle quadruples, many of which have not been determined for Apollonian quadruples, including 
the number of triangle quadruples bounded by height and by maximal element, as well as results 
about the product of elements in triangle quadruples. In Section |6} we consider the problem in 
higher dimensions. Finally, in Section|7} we collect open questions requiring additional research. 

2 Triangle Group 
2.1 Basic Definitions 

Definition 1. A triangle quadruple t = (a, b, c, d) is a quadruple of nonnegative integers satisfying 

3{a^ + b^ + c^ + d^) = {a + b + c + dY- 
Definition 2. A triangle quadruple (a, b, c, d) is primitive ifgcd{a, b, c, d) = 1. 
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/-I 2 2 2\ /I 0\ /I 0\ 

100 2 -1 22 ^_010 

010'^~00 10 '^"22-12 

VOOOl/ VOOOl/ VOOOl/ 



P' F 

f = (7,4, 3,1) f' = (7,4,9, 1) 

Figure 3: The new figure is generated by reflecting two of the segments across the side of the equilateral triangle. It is 
clear that if the original point is inside the triangle, then the new point must be outside the triangle and vice versa. 



Definition 3. The triangle group T is the subgroup ofGL{4, Z) generated by 
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2.2 General Notes About Operations 

Note that the operation generating the triangle quadruple t' = {a, b, c, a + b + c — d) from the 
triangle quadruple t = (a, b, c, d) can be characterized by t' = S4t, where t' and t are written as 
column vectors. Analogous relations hold for operations on the other elements. 

Geometrically, the operation generating (7, 4, 9, 1) from (7, 4, 3, 1) is illustrated in Figure|3} 
in the case where 7 corresponds to the square of the side length of the equilateral triangle. How- 
ever, if we remove the restrictions on order, then every quadruple actually generates four different 
configurations depending on the side length of the equilateral triangle. They can be encapsulated 
by one diagram, as illustrated in Figure |4} 



3 Reduction Theory and Root Quadruples 

Suppose we begin with the quadruple (1, 1, 3, 4). Let a reducing action be an operation that de- 
creases the sum of the elements in a triangle quadruple. We can reduce this to (1, 1,3, 1) and then 
to (1, 1, 0, 1). At this point, no further operation can reduce the sum. 

Definition 4. A triangle quadruple t = (a, b, c, d) is a root quadruple if it is not possible to apply 
a generator that reduces the sum a + b + c + d. 
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Figure 4: Segments of the same length are colored the same color. 



Lemma 1. For any triangle quadruple t = {a,b,c,d), applying the generator corresponding to 
the largest element does not increase the sum a + b + c + d. In addition, t can be reduced to 
a root quadruple t' = (0, x, x, x) (or permutations) in a finite number of operations, where x = 
gcd(a, b, c, d). 

Proof. Assume that a > b > c > d without loss of generality. From the condition, we have 

1 /3 3 

a = -[b + c + d)±d -(bc + cd + db) - -(b^ + c^ + d^). 

We first show that 



a = -(b + c + d) + ^ -(bc + cd + db) - -(b'^ + c^ + d^). 

Suppose for the sake of contradiction this is not so. If —b + c + d <0, then note that 

1 fS 3 

-(-b + c + d) < d -(bc + cd + db) - -(b^ + c^ + d^) 

1 /3 3 

^ -(b + c + d) - J -(bc + cd + db) - -(b^ + c^ + d^) < b 

=^ a<b, 

a contradiction. Therefore, —b + c + d>0. Then we have 

b^ + c^ + d^- 2bc ~ cd - 2db = b(b - c - d) + c(c - b) + d(d -b)-cd<0 

^b^ + c^ + d^- 2bc + led -2db < 6(bc + cd + db) -3 (b^ + c^ + d^) 

1 /3 3 

^ ^{-b + c + d) < d -(be + cd + db) - -(b^ + c^ + d^) 

^ a <b, 
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also a contradiction. It follows that we must have. 



1/3 3 

a = -(b + c + d) + d -{bc + cd + db) - -{b^ + c^ + d^), 

as claimed. Now, we will show that 2a>b + c + d, with equality when d = 0. Note that 



2a ^ b + c + d + ^J 6{bc + cd + db) - 3{b^ + + d^) > b + c + d, 

with equality when 

{b^ + c^ + d^)-2{bc + cd + db)=0 
^d^-2{b + c)d+{b^-2bc + c^)=0 



Assume c^0.lfd= {\/b+ ^/c)^, then clearly d> c, contradicting the minimality of d. Therefore, 
we must have d = {\/b — ^/c)^. However, then 

a=^(Z7 + c + J) = ]^{b + c + b - 2\/bc + c) = b ^- c - ^b - ^c[^ - ^c) <b, 

contradicting the maximality of a unless b = c,m which case the triangle quadruple (a, b, c, d) = 
(a, a, a, 0). The case where c = generates the same solution. □ 

In addition, it is clear from this that the nonnegativity of quadruples is preserved under this 
operation. Indeed, the product of d and a + b + c — d is 

3{a^ + b^ + c'^)-{a + b + c)^ = {a- bf + {b- c)^ + {c- af, 

the constant term of the equation 3(a^ + b^ + + d^) = {a + b + c + dY. Note that this is nonneg- 
ative, so d{a + b + c — d)>0, implying that if J > 0, then a + b + c — d>Od& well. On the other 
hand, if = 0, then a — b — c,d& proved in Lemma [T] 

Lemma 2. A triangle quadruple (a, b, c, d) can be reduced to the unique quadruple (0, x, x, x) (or 
permutations), where x = gcd(a, b, c, d). 

Proof. Suppose we begin with the triangle quadruple (a, b, c, d). Let x = gcd(a, b, c, d). Then it 
is clear that x = gcd(a, b^ c, a + b + c — d). Therefore, the greatest common divisor of the elements 
is an invariant, so by Lemma [T| by applying the generator corresponding the largest element at 
every step, it can be reduced to the unique quadruple (0, x, x, x) (or permutations). □ 

Corollary. It easily follows from Lemma^that for any two triangle quadruples t = (a, b, c, d) and 
t' = (a', b' , c', d'), if gcd{a, b, c, d) = gcd(a', b' , c', d'), then it is possible to obtain t' from t by 
applying a finite number of generators and vice versa. Therefore, all primitive triangle quadruples 
are contained in one orbit of {0, 1, 1, 1) inT. 



6 



4 Geometric Groups 



Theorem 1. The triangle group with generators S\, S2, St,, S4 is a hyperbolic Coxeter group. 
Proof. It is easy to verify that 5i, ^2, ^3, 54 satisfy 

52 =/ for {1,2, 3,4} (1) 

and 

{S.Sjf = I for / ^ J, and /, j G {1, 2, 3, 4}. (2) 
Lemma 3. The generators Si are reflections. 

Proof. The eigenvalues of 5/ are 1, 1, 1, —1. It follows that the operation corresponding to 5/ is 
the reflection with respect to the hyperplane spanning the eigenvectors v,j, vt^, vt^ of 5/. □ 

Lemma 4. For x = (a, b, c, d), each Si preserves the quadratic form 

Q{x) = 3{a^ + b^ + c^ + d^)-{a + b + c + df =xAx^, 

where 

(1 -1-1 -1\ 
-1 2 -1 -1 
-1 -1 2 -1 
V-1 -1-1 1) 

That is, Q{x) = Q{Six). 

Proof. This is easily verified by a straightforward calculation. □ 

Construct an abstract Coxeter group with A as its Cartan matrix. By Lemma [3] and Lemma |4] 
and application of basic properties of Coxeter groups, we conclude that the triangle group is the 
constructed Coxeter group. Section 5.4 in [7] proves that all relations between the matrices Si 
formally follow from ([T]) and (|2]). In other words, the geometric representation is faithful. In 
particular, since the Cartan matrix A of the triangle group has a 2 x 2 positive definite submatrix and 
negative determinant, its signature must be (3, 1), so the triangle group is indeed hyperbolic. □ 

Note that the cubic relations for Si are unique to the triangle group case, as they do not hold in 
the Apollonian circle case. The corresponding Coxeter diagram that encodes the Coxeter matrix is 
a connected graph, by ([2]). 

5 Properties of Triangle Quadruples 

5.1 The Number of Triangle Quadruples Bounded by Height 

For a triangle quadruple Q = (a, b, c, d), define H{Q) to be the height \/ a'^ + b'^ + c^ + d^ of Q. 

Theorem 2. For an integer n, define F{n) to be the number of triangle quadruples Q with height 
H{Q) < n. Then F{n) = 0{n^log^{n)). 
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Proof. Suppose that a>b>c>d without loss of generahty. Consider the substitution x = a,y = b, 
z = a + b — c,w = a + b — d. Then the quadratic form 3{a^ + b^ + c'^ + d^) — {a + b + c + d)^ can be 
rewritten as —6xy + Iz^ — 2zw + 2w^, so it suffices to show that the number of zeros of the quadratic 
form —3xy + z'^ —zw + w^ with + + z^ + < Cn^ for some constant C is 0{n^ (n) ) . It is 
clear that any solution generates a nonnegative quadruple. Indeed, if J < 0, then w > a + b = x + y, 
so z^ — zw + w^ = (z — + (3/4)w^ > {3/A){x + y)^ > 3xy, a contradiction. In addition, the 
integer property is preserved because the transformation matrix 



/I 








o\ 





1 
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-1 
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1 





-1/ 



has determinant 1 . 

Consider such a zero (w, x, y, z), so z^—zw + w^ = 3xy < Crp-, where C is a constant. Set 3xy = 
k. For a fixed k, define A (k) to be the number of integer solutions to the equation z^ — zw + w'^ = k. 
In addition, for a fixed k, it is clear that the number of integer solutions to 3xy = kis bounded above 
by d{k), the number of divisors of k. 

It follows that the number of zeros of the quadratic form —3xy + z^ —zw + w^ with x^ +y^ + 
z^ + < CrP- is bounded above by 

£ A{k)d{k). 

k<Cn^ 

Lemma 5. For any positive integer k, A{k) < 6d{k). 

Proof. Note that the function A{k) is the number of factorizations of k into a product of an integer 
and its conjugate in the ring where ^ is a cubic root of 1. This ring has 6 units (the 6th roots 
of unity). Define 5(m) : = A(m)/6. For a prime />, if p = 2 (mod 3), then p is still a prime in Z[^], 
so B{p^) = if j is odd, while B{p^) = 1 if 7 is even (unique factorization p^ = {p^ ^'^){p^ up to 
units). If p = 3 or if ;? = 1 (mod 3), then p factors uniquely in Z[^] as qq, so p^ can only factor 
as zz, where z = q\qy^\ Therefore, there are j + 1 factorizations up to units, so B{p^) = j +1. 

From these observations, it is straightforward to show that B{m) is multiplicative, implying 
that if m contains primes congruent to 2 modulo 3 raised to odd powers, then B{m) = 0, and if 
not, then B{m) = {bi + I) ■■■ {bs + I) if m is factored as p^^ ■ ■ ■ pl"''q\^ ■■■(ls\ where pi are of the 
form 3r + 2 and qi = 3 or qt = 3r + I. See [6J. In our case, it suffices to show that B{k) < d{k) for 
k = pj, a prime power. We have shown above that B{pj) =0, 1, or j + 1, and since d{pj) = j + 1, 
it follows that B{pj) < d{p^) in all cases, as desired. □ 

By Lemma [5} the number of zeros of the quadratic form —3xy + z^ —zw + w^ with x^ +y'^ + 
z^ + < Cn^ is bounded above by 

I 6(J(^))2. 

k<Crfi 

n 

Lemma 6. For any positive integer n, ^ {d{k))^ = 0(nlog^ n). 

k=\ 
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Proof. Note that 



k=l k=\i\k \<i<nk<n \<i<n^^'^y'' -1) l<Kn ^) 

j\k i<;<" 1<J<'' 1<J<« 
i\k 

d=\ " \<i'<n/d ' d=l " \<i'<n/d ' 

l<j'<n/d l<i'<n/d 
gcd(;-',/)=l 

2 



< Cnlog^n ^ "7 < Cnlog^; 



as claimed. □ 

For a generalization of this result, see Formula (1.80) in Replacing n by in Lemma |6} 
we conclude that the entire expression is bounded by 

£ 6idik)f = 0{nHog'{n)), 

k<Cn^ 

as desired. □ 

5.2 The Number of Triangle Quadruples with Specific Elements 

Suppose we want to calculate the number of triangle quadruples containing the pair of integers 
(p, q), with p,q> 0. Consider the substitution x = p, y = q, z — p + q — c, w = p + q — d, so the 
relation 3{p^ + q^ + + d^) = (p + q + c + d)^ can be rewritten as —zw + w^ = 3pq. Therefore, 
it suffices to determine the number of integer solutions to z^ —wz + z^ = 3pq- The number of 
solutions is 6B{3pq), where B{m) is defined as in the proof of Lemma [5} The closed formula 
computes the number of solutions and the solutions themselves. 

5.3 The Number of Triangle Quadruples Bounded by Maximal Element 

Theorem 3. Let max{Q) denote the maximal element in the quadruple Q. Then \{Q : max{Q) < 
n}\ = 0{nHog\n)). 

Proof. Note that max(2) < H{Q) < 2max(2), so by Theorem|2| 

\{Q : max(2) < n}| ~ |{2 : H{Q) <2n}\ = OinHog^n)) 
as well. □ 
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5.4 The Norm of Triangle Quadruples After n Operations 

Define a word w of generators Si, S2, 5*3, 54 to be a reduced word if it is not equal in T to a word 
of smaller length. Let W be the set of all reduced words. For n = Am + / with < / < 3, define the 
reduced word Rn to be Ri{S ^83528 1)'" , where Ri = I,S\, S2S1, S3S2S1 for < z < 3, respectively. 
For a matrix A, ||A||<x, is defined to be its maximum absolute row sum, so for a column vector 
V = (a, b, c, dy, ||v||<x, = max(a, b, c, d). Let r = (0, x, x, x) be any root quadruple. 

Theorem 4. For any reduced word w of length n, 

( 



\wr\ 



< ||-^n'"|U 




7 + V13 + 2 




Proof. Write w = Si^Si^_^ ■ ■ ■ 5,j . Let w" = wr and r'^ = Rnr. Suppose that the elements of w" and 
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►•4 '►'1 ^ '4 ~ ' 1 (4) 

by induction on n. The base case with n = 1 is trivial. For the inductive hypothesis, assume that 
the claims hold for n = k—\. To prove that ([3]) holds for n = k, note that 
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To prove that Q holds for n = k, consider two cases. If w\ 
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- ^« _ 

~ ' 4 ' 1 • 

This completes our induction. Note that the proof is easily derived from the proof of Theorem 7.1 
in Q with a few modifications. 

To bound | |7?„r| |oo, consider the characteristic polynomial of 5453S251 , l—lt — I5t^ — + 1^. 
Note that the growth of (545352S1)"' is bounded by y"\ where 



is the largest eigenvalue of 54S352S1 and the largest root of its characteristic polynomial. □ 



5.5 The Number of Triangle Quadruples Generated After n Operations 

Theorems. Let Wn be the set of reduced words with length not greater than n. Then \Wn\ = 0{X"), 
where A = ^'^'f^ . 

Proof By [|T1, the growth coefficients G„ of the triangle group, computing the number of reduced 
words with length n generating distinct matrices, are generated by the recurrence Go = 1, Gi = 4, 
G2 = 12, and G„ = 2G„_i + 2G„_2 — 3G„_3 for n > 3. The solution for this recurrence is 

Therefore, if we let A denote then G„ ~ A". It follows that 

= Go + Gi+G2 + --- + G„^ 1+A+A^ + ---+A"= "^.^ ~^ =0(A"), 

A — 1 

as desired. □ 

Theorem 6. For any root quadruple r of the form (0, x, x, x), ci?i"/n^ < \Wnr\ < C2?i''for some 
constants ci and C2- That is, | W„r| and \W„\ differ asymptotically by a factor that is 0{n^). 

Proof. The upper bound is trivial by Theorem [5} To prove the lower bound, note that the stabilizer 
of r in r is Tq := S3 x I? (where S3 denotes the symmetric group in three items), the subgroup 
generated by 52, ^3, 54. Therefore, if w and w' G Wn are in the same coset of Tq, then w' = wu for 
some u E Tq. It follows that u = w^^w', so the length of u is at most 2n. But by fTl, the growth 
series of Tq, which has three generators, is tj^^^ from which it can be computed that the growth 
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coefficients G„ are generated by Go = 1 and G„ = 3n for n> 0. Therefore, the number of elements 
in Tq of length at most 2n is 

1+3(1 + 2H \-2n) =6n^ + 3n+l <kn^, 

for some constant k. For any c G W„r, define N{c) := \{w E Wn :wr = c}\. Then 

\Wn\= £ Nic) < \W„r\kn^, 

ceW„r 

SO I W„r| > \Wn\/krP-, which implies the lower bound. 

□ 

5.6 Elements in Triangle Quadruples 

Let 0(4) := 0(4, A, C), the group of matrices that preserve the quadratic form A. Similarly, let 
50(4) :=50(4,A, C). 

Theorem 7. For every triangle quadruple Q = (a, b, c, d), define Cc{Q) to be the number of prime 
factors (with multiplicities) in a-b ■ c ■ d. There are infinitely many Q such that Oc{Q) < kfor some 
constant k. 

Proof. We will invoke Theorem 3 from [l4l|, which states that the result is true if the Zariski closure 
r of r in the complex algebraic group GL(4, C) is a semisimple algebraic group. By this theorem, 
it suffices to verify that T is semisimple. 

Lemma 7. The Zariski closure of the triangle group, T, is 0(4). 

Proof. Consider the subgroup T123 = S3 x (where S3 denotes the symmetric group in three 
items) generated by 5i, ^2, ^3, containing the translation matrix 
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It is straightforward to verify that 
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It follows that Bi E Lie(r), along with the conjugates of 5i by T and their brackets. It is 
straightforward to verify that the six matrices 



-1 2 -1 -1\ /2 -1 -1 2 

r- D c-1 _ / -1 2 -1 2 \ c R c-1 _ / 2 -1 -1 -1 

0/ Vo 



{SiBiS^')iS4B,S^') - iS4BiS^')iSiB,S;'] 



3 -6 12 -6 
12 3 -6 -6 
-6 12 3 -6 

-9, 



3 12 -6 -6 

(525i52 ^)(545i54 1) - (S45i54 l)(525i52 ^) = | -6 ^3 -t 

-9, 

(30 12 12 -24' 
12 -6 -6 -6 
36 -18 -18 -If 



[5254^1 (5254)"^] (54^154 1) - (545154^) [5254^1 (5254)" 



-6 12 -6 -6 
12 30 12 -24 
-6 12 -6 -6 
-18 36 -18 -18, 



are linearly independent. These six linearly independent elements of Lie(r) determine a basis, im- 
plying that they span Lie(50(4)) (with dimension 6), which equals Lie(0(4)), so Lie(0(4)) C 
Lie(r). In addition, by Lemma Q T C 0(4), so Lie(r) C Lie(0(4)). Therefore, Lie(r) = 
Lie(0(4)). This implies that T is either 50(4) or 0(4). However, T cannot be 50(4) because 
T contains matrices of determinant —1 (5i, for example), so T = 0(4), as desired. □ 

It is well-known that 0(4) is semisimple, and this completes our proof. □ 



6 Extension to Higher Dimensions 

Lemma (a) in [5 J proves the following result, which extends the equation for triangle quadruples 
to higher dimensions. Here, we present an alternative proof. 

Theorem 8. For a simplex T in n dimensions and an arbitrary point P, if we let ao, ai, . . . , an+i, 
denote the square of the side length ofT and the squares of the distances from P to the vertices of 
T, then ao, ai, satisfy (n-F l){a}^^-a\-\ ha^^J = (ao + ai H hQw+i)^. 

Proof. Consider a point P inside a simplex in n dimensions, configured so that P is at the origin. 
Then the vertices of the tetrahedron are the vectors vi, . . . , v^+i. Set = k; — and a,- = |v,p, 
so 2(v/, Vj) = aj + aj — a^. Since v,- G R" and there are n -t- 1 of them, they are linearly dependent, 
which implies that the Gram determinant (the determinant of the matrix with entries (v,, vj)) is 
zero. Therefore, det(G) = 0, where Gij = ai + aj — ao, for i, j = I, n + I, i ^ j, and Gu = 2a,-. 
Let a be the vector with entries a,-, for z = 1, . . . , n + 1, and let u be the vector whose entries all 
are 1. Then, for any column vector x = {xi, . . . , Xn+i) E R"+\ we have (G — aoI)x = (u, x)a + 
(a — aou, x)u. Since the image of G — aol has dimension 2, it follows that G — aol has rank 2 and 
consequently, an n — 1-dimensional kernel K. Set ei = a, ^2 = u, and select a basis ^3, . . . , Cn+i of 
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K. In this basis, all the rows of the matrix G — oqI are zero, except the first two, and the upper left 
2x2 block contains the elements 

(u, a) (u, u) 

(a — aou, a) (a — aou, u). 

The characteristic polynomial of G — qqI is thus the characteristic polynomial of this 2x2 matrix 
multiplied by or 

(z^ - (2a - uqu, u)z + (u, a) (a - uqu, u) - (u, u) (a - uqu, a) ^ . 

The characteristic polynomial of G is obtained from this by replacing z by z — ao, or 

((z-^o)^- (2a-(2ou, u)z+ (u, a)(a-aou, u) - (u, u)(a-aou, a))(z-«o)"~'^- 

Therefore, substituting z = 0, we calculate 

det(G) = (-l)«+Met(-G) 

= (-1)"+^ (-flo)"~^ («o + <^o(2a - aou, u) + (u, a) (a - aou, u) - (u, u) (a - aou, a)) 
= aQ~\ao + ao(2a — aou, u) + (u, a)(a — aou, u) — (u, u) (a — aou, a)) 

= aQ~^[{ao + ai-\ ha„+i)^ - (n + l)(ao + af H 

Since det(G) = 0, we conclude 

(flO + «i H 1- «n+i )^ - (« + 1 ) (oq + «i H 1- ) = 0, 

as desired. □ 
Note that the general relation simplifies to 

n n / " \ ^ 

na^+l-2a„+i J2 + +1)1^ =0' 

(=0 (=0 \i=0 / 

so we can reduce as follows 

2 " 

(ao, ai, Un+i) {ao, ai, - V a, — a„+i). 

".to 

Notice that in dimensions n > 2, the analogous operations do not preserve the integer property of 
the elements. 

7 Open Questions 

1. Beginning with a specific root quadruple, is it possible to calculate the asymptotics of the aver- 
age value of the maximum element in the triangle quadruple obtained after n operations as n goes 
to infinity? 
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This problem follows naturally from our discussions of the number of triangle quadruples char- 
acterized by maximal elements and by height, as well as the growth rate of the triangle group. 
Intuitively, the average value of the maximal element after n operations should grow exponentially. 

2. Given a triangle quadruple, is it possible to calculate how many reduced words of minimal 
possible length can generate it? 

This is related to the previous question in that it requires an analysis of the properties of trian- 
gle quadruples after n operations. 
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